Anderson universality in a model of disordered phonons 
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We consider the localization properties of a lattice of coupled masses and springs with random 
mass and spring constant values. We establish the full phase diagrams of the system for pure 
mass and pure spring disorder. The phase diagrams exhibit regions of stable as well as unstable 
wave modes. The latter are of interest for the instantaneous-normal-mode spectra of liquids and 
the nascent field of acoustic metamaterials. We show the existence of delocalization-localization 
transitions throughout the phase diagram and establish, by high-precision numerical studies, that 
the universality of these transitions is of the Anderson type. 

PACS numbers: 63.20.D-, 63.20.Pw, 63.50.-x 



Coherent wave phenomena in disordered systems are a 
recurring theme of modern physics. For condensed mat- 
ter at the quantum scale, Anderson localization [l|, |2| 
has recently been re-modernised by a series of beauti- 
ful experiments for Bose-Einstein condensates and the 
exponential decay of the localized waves has been di- 
rectly measured 0] . Results agree very well with recent 
spatially- resolved studies in semiconductor systems [jj. 
Similarly, light localization continues to be at the fore- 
front of many research activities [5|. For classical waves, 
localization phenomena have an equally impressive his- 
tory [6| and recent ultra-sound propagation experiments 
[7| can now probe the spatial and multifractal structure 
of states close to the Anderson transition. Phonon local- 
ization, i.e. absence of diffusion of acoustic or vibrational 
degrees of freedom has been addressed first in the seminal 
paper by John et al. [8| and thereafter in connection with 
the low-temperature thermal properties of glasses and the 
related enhancement of the vibrational density of states, 
the so-called "boson peak" J9l4l3j. Also the localization 
properties of the Hessian Matrix of the potential-energy 
landscape of liquids and glasses — "instantaneous nor- 
mal modes" — have been investigated and related to the 
liquid-glass transition [14|. 

A very interesting new avenue of research has opened 
up recently due to the realisation that the seminal work 
on electromagnetic metamaterials [l5| has a companion 
in acoustic systems as well |16j . Hence hitherto unex- 
plored and deemed unphysical regions of the phase dia- 
gram for disordered vibrations — those with apparently 
negative masses and spring constants — are now recog- 
nised to be of considerable interest for metamaterial ap- 
plications and offer an entirely novel perspective of An- 
derson localization. 

In this paper, we address the localization properties of 
a simple cubic lattice of particles with varying mass nii 
and varying nearest-neighbour harmonic force constants 
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We present for the first time the complete phonon 



stable regime. We achieve this by using high-precision 
transfer- matrix methods (TMM). 

For simplicity we deal with "scalar displacements" 
Ui{t) ~ Ui{uj)e^'^^ [11| obeying equations of motions in 
the frequency domain 
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where the sum runs over the 6 nearest neighbours of a 
lattice site i. 

In our calculations we use two types of disorder: In 
the case of mass disorder we allow the masses rui to be 
uniformly distributed in the interval [r?T — Ato/2, m -I- 
A?7i/2] with TO = 1 and keep kij ~ k = 1 constant. In the 
case of spring disorder we keep m^ = to = 1 constant and 
distribute the kij uniformly in the interval [1 — A/c/2, 1 -|- 
Afc/2]. If the width of the distributions Am and Afc 
exceeds 2, negative masses/spring constants appear and 
mimic the unstable part of the Hessian of a liquid or the 
harmonic properties of an acoustic meta-material. 

The present model can be related to Anderson's [l| 
tight-binding model for electron localization with Hamil- 
tonian H = J2i l*)^j(*l + Sii K)*JiOI- For mass disorder 
we set tij = kij = const. = 1 and obtain the relations (E 
is the quantum energy) 



E = -LU^ 
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e,; = w^ (mi - 1) , 



(2) 



while for spring disorder (rUi = m = 1), we have 



E: 
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we can reuse many 
ii. In 
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localization diagram for such a system, including the un- 



With the relations ^ and Q 
of the results for the Anderson model 
particular, the analogy establishes the existence of 
localization-delocalization transitions for the vibrational 
mass-disorder model. We introduce the electronic width 
parameter W = AttiJo;^! for the width of the electronic 
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FIG. 1: (Color online) Phase diagrams of (a) mass disorder with disorder parameter Am and (b) spring disorder with disorder 
parameter Afc versus squared frequency oj^. Shaded lines in (a) denote the critical region obtained from the transformation of 
the electronic phase diagram of [l7|] via Eqs. ([2]). Light shadings denote localized regions, dark shadings regions beyond the 
band edges. Open blue diamonds denote transition point determined by finite-size scaling. Solid red circles denote estimated 
transition points (see text) and green crosses denote maxima in the DOS, divided by lo'^ ("boson peaks"). The large open circles 
denote the locations of the critical eigenstates in Fig. [21 The horizontal dash-dotted lines indicates the border between stable 
and unstable regions, the dotted lines denote to^ = and 6. The red line for (b) marks a region of low numerical accuracy 
where the phase boundary is less well determined. Error bars for TMM data as described in the text. For the boson peak, 
errors denote width of the peak at 95% of its height. Insets: blow-ups of the phase boundaries in the stable regions. 



disorder distribution. This allows us to estimate the vi- 
brational mass disorder phase boundary from the elec- 
tronic on-site disorder phase boundary [17| . 

In Fig. [Ija) we show the estimated mobility edges for 
the case of vibrational mass disorder. The phase diagram 
is intriguing in many respects. We first note that the re- 
gion for w^ > 6 corresponds to the E < region in the 
Anderson model and similarly w^ < 6 is associated with 
E > 0. The much studied centre of the band at £■ = 
for the Anderson model becomes the much less distinct 
u;^ = 6 line in Fig. [TJa). For w^ < 6, we see that much 
of the extended phase belongs to the region of possible 
negative masses with Am > 2. Hence this region cor- 
responds to a well studied counterpart in the Anderson 
case. Furthermore, the E > region also extends into 
negative values of w^ and gets transformed into a much 
reshaped form for o;^ < 0. The particular form of this 
puddle of extended states, towards the cj^ = axis, is 
driven by the so-called re-entrant behaviour for the An- 
derson model 17l.ll9l|. Similarly, the re-entrant behaviour 



at w > 12 can be traced to the corresponding re-entrant 
shape of the mobility edge a,t E < —6. As we will show, 
these extraordinary mobility edges and hence the phase 
diagram for the mass disorder case are indeed confirmed 
by our direct high-precision numerics. 

For spring disorder, ([3]) corresponds to a disorder dis- 
tribution consisting of the sum of 6 independently chosen 
random numbers. Even when each kij is chosen accord- 
ing to the uniform distribution as above, the resulting 
distribution of e^ has not previously been studied for the 



Anderson model — although studies with pure hopping 
disorder exist [20[ — such that there are no phase dia- 
grams to compare with. We have determined this phase 
diagram using our high-precision TMM. 

In the TMM a quasi-onedimensional bar with fixed 
cross section M x M for lengths L ^ M is considered. 
Equation ([1]) is then rearranged into a form in which the 
amplitudes of vibration for a cross-sectional sheet can 
be calculated solely from those in the preceding cross- 
sectional sheet. If we denote the i-th M x M sheet of 
displacements by C/,, Eq. ([T|) can be expressed recursively 
as {Ui+i,Ui) = Ti-{Ui,Ui-i) with suitably defined trans- 
fer matrix T^ ^ . The standard TMM for the Anderson 
problem [21| is then used to calculate the Lyapunov expo- 
nent Am of the mapping via T^'s and the reduced decay 
length Am = Xm/M- 

We performed TMM calculations at various values of 
Am, Afc indicated in the phase diagrams in Fig. [T] For 
every disorder value, we calculate Am for a range of fre- 
quencies and system widths M = 6,8, 10 and 12 to a rel- 
ative error of 0.1%. The transition is initially estimated 
as the frequency where the M = 12 and M = 10 lines 
cross. The error shown is the difference with respect to 
the frequencies where the M = 12 and the M = 6 data 
cross. We see from Fig. [ija) that these rough estimates 
are in excellent agreement with the phase boundaries as 
established from the analogy with the Anderson transi- 
tion. In Fig. [TJb) the same method is used to obtain the 
phase diagram for the spring disorder. In addition, we 
performed studies at larger system sizes, up to AI — 20 



with 0.1% error, for the delocahzation-locahzation tran- 
sitions at three representative regions in the phase dia- 
grams indicated by open diamonds in Fig. [T] We found 
clear transitions from extended behaviour, with increas- 
ing Am values for increasing M, to localized behaviour, 
where Am decreases when M increases. The transition 
frequencies obtained are in excellent agreement with the 
estimates by the method described above. 

For spring disorder we see that in the "central region" 
around lu^ = 6 states remain extended up to the largest 
considered spring disorder Afc = 10. This is similar 
to the electronic case with pure hopping disorder [20| 
where even very strong hopping disorder does not lead 
to complete localization close to -E = 0. Whether the 
re-entrant behaviour of the phase border above the in- 
stability line for w^ > 8 is genuine, remains to be de- 
termined by higher-precision calculation. We note that 
it coincides with the vanishing of the boson peak. Such 
a re-entrant regime is absent in the off-diagonal Ander- 
son system. For the spring-disorder model it could sig- 
nify a combined effect of localization and instability. For 
u!^ < (and Afc > 4), we observe an even larger area of 
extended states than for mass disorder. The localization- 
delocalization transition on the co^ < side is very similar 
to that observed in the instantaneous-normal mode spec- 
tra. The significance of the delocalizcd unstable modes to 
the energy landscape of a liquid remains to be discussed. 

We find that both for mass and spring disorder, the 
w = hydrodynamic mode remains extended regardless 
of the disorder strength. This is in agreement with pre- 
vious studies in one- and two-dimensional systems |13l |. 
What is also common to both mass and spring disorder 
is the observation of very strong shifts of the crossing 
points of Ajvf when changing M. Such a behaviour is 
to be expected, however, since we are effectively dealing 
with transition regions in the vicinity of the band edges 
where density-of-states effects can dominate the scaling. 
This is again similar to the situation for the electronic 
case where the transition at the mobility edges ior E ^ 
is also more difficult to study [2^, \2M ■ 

We have also computed the density of states g{uj) of 
our models by directly diagonalizing the dynamical ma- 
trix. We have divided g{uj) by w^ in order to detect 
maxima which correspond to "boson peaks" , which have 
been discussed in the literature. In the stable regime be- 
low the boson maxima the vibrational excitations are es- 
sentially wave-like excitations. For mass disorder we find 
such maxima within the whole phase diagram between 
a; = and the mobility edge as shown in Fig. [Ija). For 
spring disorder the boson peaks disappear slightly above 
the instability line Afc = 2 (cp. Fig.[ljb)). We have some 
evidence from our density-of-states calculations and from 
the form of the wave functions that indeed wave-like exci- 
tations persist in the instable region of the mass-disorder 
model, whereas there is no evidence for wave-like exci- 
tations in the spring disorder model for Afc > 2.5. This 



Am 


M 


OJ^ 


u^i 


V r. 


1.2 


8 to 20 


12.15 to 13.1 


12.68 ± 0.06 


1.57 ±0.14 0.84 


4 


8 to 20 


3.75 to 4.25 


4.13 ±0.04 


1.57 ±0.08 0.99 


9 


8 to 20 


-1.65 to -1.5 


-1.62 ±0.04 


1.57 ±0.41 0.87 


Afc 


M 


C.2 


-c^ 


V Tq 


1 


10 to 20 


12.48 to 12.6 


12.527 ± 0.003 


1.58 ±0.05 0.62 


10 


6 to 16 


18.8 to 20.3 


19.75 ±0.05 


1.51 ±0.08 0.84 


7 


8 to 20 


-3.5 to -2.75 


-3.33 ±0.11 


1.59 ±0.29 0.51 



TABLE I: Values of critical parameter cu^ and i/ for spring 
and mass disorder computed from FSS performed in the given 
M and uj^ ranges. The goodness-of-fit parameter Vq is also 
shown for each fit. Errors denote the standard error-of-mean 
estimates. 



striking difference may be due to the fact that for mass 
disorder the disorder fluctuations are suppressed by a 



factor 



leaving a slightly disturbed spectrum of the 



simple cubic lattice in the small-w regime. For spring dis- 
order the band character is obviously destroyed already 
for small values of Afc. 

We turn now to a discussion of our high-precision de- 
termination of the critical parameters. In order to as- 
certain the existence of a divergent correlation length 
^(w^) oc \uj^ ~ uj1\~" at ujI with critical exponent i/, we 
need to proceed, as usual in the electronic case, via a 
finite-size scaling (FSS) procedure [2J]. The FSS includes 
corrections to scaling which (i) account for the nonlincar- 
ities of the Am, Afc dependence of the scaling variables 
(relevant scaling) and (ii) for the shift of the point at 
which the Aa/(w^) curves cross (irrelevant scaling). This 
analysis is by now standard and we refer to the litera- 
ture for details of when fits are acceptable as stable and 
robust as well as for error estimates via Monte-Carlo ap- 
proaches [2J, |25[ . Details for the chosen expansions in the 
present case can be found in Ref. 118[. In Tab. HI we show 
the results for the high-precision FSS analysis at the 6 
representative disorders. We find that in all cases, a con- 
sistent, robust and stable fit can be found with quality- 
of-fit parameter Tq larger than 0.1. As our results show, 
we find that the critical exponents for both mass and 
spring disorder in the stable, unstable and — w^ regions 
of the phase diagram agree with each other within the 
error estimates. In addition, they agree equally well with 
current estimates of the correspon ding ex ponent for the 
Anderson model of localization 122. 24- 



261. Therefore we 



conclude that the scalar model of lattice vibrations stud- 
ied here falls into the universality class of the Anderson 
transition 127 



28| 



In Fig. [2] we have represented critical amplitude distri- 
butions of the states marked with black circles in Fig. [T] 
We see that the structure for mass disorder is quite dif- 
ferent from that of the spring disorder. An evaluation of 
the multifractal properties of these states might probably 
exhibit different spectra for the two models. 




FIG. 2: (Color online) Schematic representation of critical 
amplitude distributions obtained from exact diagonalization 
for (a) Am = 4 at o;^ = 4.134 and (b) Ait = 1 at cj^ = 12.526. 



All sites with 



E=. 



> 1 are shown as small cubes and 



those with black edges have ^ 



> vTOOO. The color 



scale distinguishes between different slices of the system along 
the axis into the page. 



In conclusion, our results show that a disordered scalar 
phonon model exhibits all the rich features of the Ander- 
son localization-delocalization transition. While the crit- 
ical exponents are universal and of Anderson type the 
mass-disorder and spring-disorder models exhibit com- 
pletely different localization phase diagrams. The re- 
entrant behaviour of the mass-disorder system — inher- 
ited from the Anderson model with on-site disorder — 
is present both on the stable and unstable side of the 
phase diagram. The spring-disorder phase diagram is 
dominated by delocalized states. Localized states exist 
on both sides near the band edges. In the stable regime 
Am < 2, Afc < 2 localized states exist only at the upper 
band edge in agreement with earlier investigations |8l.lll 
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